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We analyze both the s- and p-wave collision induced frequency shifts and propose a over-7r pulse 
scheme to cancel the shifts in optical lattice clocks interrogated by a Rabi pulse. The collisional 
frequency shifts are analytically solved as a function of the pulse area and the inhomogeneity of 
the Rabi frequencies. Experimentally measured collisional frequency shifts in an Yb optical lattice 
clock are in good agreement with the analytical calculations. Based on our analysis, the over-7r 
pulse combined with a small inhomogeneity below 0.1 allows a fractional uncertainty on a level of 
10 -18 in both Sr and Yb optical lattice clocks by canceling the collisional frequency shift. 
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Ultracold atoms trapped in an optical lattice have 
served as ample frameworks for high-precision atomic 
clocks 5H1, quantum information processing [1, @, 
and for quantum simulators with Bose-Einstein conden¬ 
sates [71, Isj. Many atoms in the same optical lattice site 
inevitably lead to collisions among atoms; therefore, ef¬ 
forts to understand and control the collision dynamics are 
crucial to improve performance levels. For optical lattice 
clocks, the simultaneous interrogation of many atoms is 
a key feature in contrast to ion clocks, and a reduction of 
the collision shift is one of the most important factors 
to reach relative frequency uncertainty level of 10 -18 . 
To overcome collisional shifts, the suppression of these 
collisions by strong atom-atom interactions in a Sr 2D 
optical lattice [12J and precise measurements of the colli¬ 
sion shift per atom number in a Sr ID optical lattice with 
Rabi spectroscopy 0 were reported. Collision-shift can¬ 
cellation has also been demonstrated with Ramsey spec¬ 
troscopy 13, 15] and with Rabi spectroscopy [9 . fill. Tiul] . 
Recently it has been shown that the inhomogeneity of 
Rabi frequencies allows s-wave collisions j9H12l and that 
p-wave collisions are even dominant in Yb and Sr lattice 
clocks [ii Izl. However, in Rabi spectroscopy, the can¬ 
cellation of the shift were investigated by only terms of 
s-wave collisions in the case that the initial state is the 
excited state. The p-wave collisions gives the dramatic 
difference of the shift with respect to the s-wave colli¬ 
sion case, especially when the initial state is the ground 
state. Nowadays the Rabi line shape distortion in the p- 
wave dominant case has been studied in aspect of observ¬ 
ing many-body interactions 0- In contrast to Ramsey 
spectroscopy, the cancellation of the p-wave collision in¬ 
duced shift, which is important in optical lattice clocks, 
has not been explored and analytically analyzed in Rabi 
spectroscopy. 

In this Letter, we analyze both s- and p-wave cold col¬ 
lisions and propose a over-7r pulse interrogation scheme 
for canceling the collision shift in optical lattice clocks 
using Rabi spectroscopy. We apply our analysis to the 
experimentally measured collisional frequency shift in 


an Yb optical lattice clock. Although shift and uncer¬ 
tainty on the 10 -17 level are reported due to operation 
at slightly off of collision-shift-free conditions, it is shown 
that shift and uncertainty on the 10“ 18 level can be 
reached under experimentally achievable conditions, such 
as a small inhomogeneity below 0.1 in the Rabi frequen¬ 
cies of trapped atoms and a small degree of lock servo- 
fluctuation (ztl %), for both Yb and Sr lattice clocks 
when using an optimum over-7r pulse. 

In actual ID optical lattice cases, more than two atoms 
exist in a disk. Although a many-body Hamiltonian 
is required to describe the ID case Gim , the sim¬ 
ple two-atom four-level Hamiltonian shown in Eq. m 
with effective collision-interaction frequencies can pre¬ 
dict many-atom cases in the weak atom-atom interac¬ 
tion regime by adopting a mean field picture [13]. In 
order to calculate the collisional frequency shifts in Rabi 
spectroscopy, we solve the Schrodinger equation with the 
four-level Hamiltonian on a two-atom basis, in this case 
singlet(| eg~) = ^ [| eg ) _ \ ge )]) a nd triplet states^#), 

|ee), and | eg + ) = [\eg] + \ge)\), 
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where V gg and V ee , 

and V eg and U eg 

correspond to the 


frequencies characterizing the p-wave collision shifts of 
the triplet states and the s-wave collision shifts of the 
singlet state, respectively. These frequency shifts are de¬ 
termined by the scattering lengths and the atomic den¬ 
sity. 5 denotes the laser detuning with respect to the 
bare atom resonance frequency, oo eg . is the frequency 
error of the clock laser. Cl is the thermally averaged Rabi 
frequency of atoms in vibrational states in an optical lat¬ 
tice and A Cl is the RMS spread in the Rabi frequency, 

A Cl — yj (ClV) T — Cl 2 , where is the thermal aver¬ 

age of the square of the Rabi frequencies. We define the 
inhomogeneity of the Rabi frequency as 7 = A Cl /Cl. 
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FIG. 1: (Color online) (a) (Left) Energy level diagram of 
ultracold two-atom collisions on triplet and singlet bases. 
(Right) Energy level diagram on a single-atom basis, (b) 
Excited population as a function of the laser frequency de¬ 
tuning and pulse area, obtained by projecting the singlet and 
triplet populations onto the single-atom basis. The white 
solid line depicts the peak of the transition spectrum, (c) 
Time evolution of the probabilities of singlet and triplet 
states. The red shaded (or hatched) region represents the 
positive (or negative) time-averaged population difference, 
fo N s(.t) - N e (t)dt/T > 0 (or / Q T N a (t) - N e (t)dt/T < 0). 


The populations of the triplet and singlet states are 
projected to a single-atom basis, \e) and | g), as shown 
in Fig. U-(a), to obtain the e-g transition spectrum 
which is the expectation value of the projection opera¬ 
tor, P e = \ee) (ee| + \ \eg + ) (eg + 1. The collisional fre¬ 
quency shift of the clock transition is determined by 
measuring the center of the e-g transition spectrum, as 
shown in Fig. [l]-(b). As intuitively expected and shown 
in more detail below, the collisional frequency shift de¬ 
pends mainly on the time-averaged population difference 
(i.e., / 0 T N g (t) — N e (t)dt/T, where T is the pulse dura¬ 
tion of a Rabi pulse). This is shown in Figs. [I]- (b) and (c), 
where the X point (or Y point) with a positive (or neg¬ 
ative) shift shows a positive (or negative) time-averaged 
population difference. 

With the help of the second-order perturbation theory 
under the condition of Q Vij =ej9 , U eg , S, we can obtain 
the analytical formula of the collisional frequency shift as 
a function of the pulse area A = QT. The collisional fre¬ 
quency shift is represented by the sum of a homogeneous 
p-wave, an inhomogeneous p-wave, and the s-wave colli- 
sional frequency shifts, 6 coi = S^ ol + 6$\ ol + 5^ col . 

The homogeneous p-wave collisional frequency 


shift (HP CFS) is the dominant term among the three 
types of collisional frequency shift. It comes from 
the g-e population difference due to the time-varying 
probabilities of the triplet states assisted by Cl. It can 
be written as 

4cL * -^Ve 5 - V gg )a{A) + (Vee - V gg )f3(A), (2) 

where a(A) = (-1 + s^) /(2£(A)), p(A) = 

[ 2s A ( T 2) + - 3)]/£(A), and £pL) = 

^sm_^/2) _ Because the time-averaged population 

difference is given by Jq N g (t) — N e (t)dt/T ~ sin A/A 
in the perturbative regime, it is revealed that the ho¬ 
mogeneous p-wave collisional frequency shift depends on 
the time averaged population difference, S^ ol ~ C\ + 
C 2 f 0 T N g (t) — N e (t)dt/T, as shown in Fig. [2]- (a). 

The inhomogeneous p-wave collisional frequency 
shift ( IPCFS) is determined by the g-e population differ¬ 
ence resulting from the leakage of the triplet-state popu¬ 
lation to the singlet state due to AQ. It can be written 
as 

4hloi « 1 2 m)(V eg - V gg ) + v(A)(V ee - V gg ) (3) 

-x(A)V gg ], 

where 6(A) = [—3 + 4 A 2 + 2A 2 cosA + 3cos2A — 
A 3 sin A] / (4A 2 £(A)), r](A) = [75 - 6A 2 + 24(-3 + 
A 2 ) cos A —3 cos(2A)— 60Asin A+7A 3 sin A]/(24A 2 ^(A)), 
and x(A) = [cos 2 f - (^ + ^)^]/^(A). The IP CFS 
increases with the pulse area, dominated by 6(A). 

The singlet-state | eg~) populated according to the 
nonzero AQ undergoes an s-wave collision. The s-wave 
collisional frequency shift (SCFS) is written as 

7 2 x(A)U eg . (4) 

The SCFS is 0 at A = 7r, which indicates that the time- 
averaged population difference is 0 H10- This is the 
smallest term among the three types of collision shifts. 

We normalize the collisional frequency shift to 5 co i = 
(V 6c Cv ) verify the zero-crossing behavior of the col¬ 
lisional frequency shift independent of V eg — V gg . Each 
normalized collision shift is plotted in Fig. [2]- (a). The 
negative HP CFS in the over-7r pulse area regime is the 
heart of the cancellation of the collisional frequency shift 
because the IP CFS is always positive (see Fig. [2]- (a)). 
To realize the cancellation of the negative HP CFS and 
the positive IP CFS, a sufficiently small inhomogeneity is 
required. 

The HP CFS can be positive in the entire range of the 
pulse area in the case of a large V ee — V gg values; hence 
the cancellation of the collisional frequency shifts is not 
possible in such a case. To check the necessary condition 
of V ee — V gg for cancellation, i.e., S^ ol < 0, the HP CFS is 
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FIG. 2 : (Color online) (a) Normalized HPCFS and IPCFS as 
a function of the pulse area. The blue solid line and the black 
dashed line depict the HPCFS , S^ ol , w ^en a — 0.4 and the 
time-averaged population difference, respectivley. The red 
dashed-dotted line and the purple dotted line represent the 
IPCFS , b[hcoV when 7 = 0.1 and 7 = 0.3, respectively, (b) 
Normalized collisional frequency shift as a function of a and 
the pulse area. 


(a) Sr like, g=0.4 
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FIG. 3: (Color online) (a) and (b) Normalized collisional fre¬ 
quency shift S co i as a function of the pulse area and the inho¬ 
mogeneity in Sr and Yb. (c) and (d) Normalized uncertainty 
u co i along with the zero shift curve in Sr and Yb. It is as¬ 
sumed that SA = 0.1 % and £7 = 5 %. 


calculated as a function of the pulse area and the ratio be¬ 
tween the collision interaction frequencies a = Yl e ~^ 99 , 

v eg Vgg 

as shown in Fig. [2]-(b). The small value of a < 0.63 is a 
necessary condition for the cancellation of the collisional 
frequency shift. Fortunately, the cancellation of the colli¬ 
sional frequency shift is possible in both the Sr (a ~ 0.4) 
and Yb (a ~ 0.1) optical lattice clocks 13l fl 8 |. 

In the beyond-perturbation regime, i.e., Cl < Vij= e , g , 
the Rabi spectrum is distorted and the shift starts to 
stray from the mean-field approach jl 8 :]. The shift is not 
clearly defined in the distorted Rabi line shape. In an 
atomic clock based on the Rabi spectroscopy, its opera¬ 
tion in the perturbative regime is appropriate. 

We calculate the total collisional frequency shift as a 
function of the pulse area and the inhomogeneity in both 


Yb and Sr atoms using Eqs. ([2])-(j4]). Figures [3]- (a) and 
(b) show the existence of the zero-collisional frequency 
shift in the over-7r pulse regime. 

Figures [3]- (c) and (d) show the normalized uncertain¬ 
ties in the S co i(A , 7 ) = 0 curves of Sr and Yb. Fluc¬ 
tuations in the atomic density, pulse area, and inho¬ 
mogeneity are the main origins of the residual colli¬ 
sional frequency shift change. The normalized deviations 
of the residual collisional frequency shift changes (i.e. 
normalized uncertainties) are calculated by u co i = 

\J (<L ;^) 2 + 0. In the col- 

lision shift-free 5 co i(A , 7 ) = 0 curve, the uncertainty 
stemming from the fluctuation of the atomic density, 
S co i^-, disappears. Table U shows the conditions for 
u co i < 2 x 10 -3 , as depicted by the O-A line in Fig. [3j 
V eg — V gg under the Hz level can provide 5 x 10“ 18 uncer¬ 
tainty in this condition. The atomic density should be 
on the order of 10 10 /cm 3 ; this can be obtained by the 
large volume optical cavity of an optical lattice without 
sacrificing the S/N ratio. 


TABLE I: Collision-shift cancellation condition for u co i < 2 x 
10 “ 3 . 


Atom 

a A 

7 |v;,, 

— Vgg | for a 5 x 10 18 




(Hz) 

Sr 

0.4 A < 1.257r 

7 < 0.09 

< 1.07 

Yb 

0.1 A < 1.05tt 

7 < 0.11 

< 1.30 


To apply our analysis to experiments, we measure the 
collisional frequency shifts of an m Yb optical lattice 
clock as a function of the pulse area in two inhomo¬ 
geneity cases, 7 = 0.23 and 7 = 0.15, the latter be¬ 
ing the lowest in our current experimental setup. The 
details of our Yb optical lattice clock are described in 
the literature |20|, |2l|. The potential depth of the ver¬ 
tically oriented ID optical lattice is Uo = 185i% and 
the trap frequencies are v z ~ 55 kHz and v r ~ 202 
Hz, as measured by the sideband spectra [22]. Nearly 
50000 atoms are trapped in the optical lattice. The 
atom density here is roughly po = 4.3 x 10 10 atoms/cm 3 
where Vij= e , g and U eg are sub-Hz levels. To measure 
the collisional frequency shift at the given atom den¬ 
sity, the interlacing between the high (No) and low atom 
numbers(0.577Vo) is repeated. The final excited fraction 
is fixed at 90% (± 6 %) of the maximum excited fraction by 
feedback to the computer-controlled AOM for the clock 
laser (a clock laser lock servo) [23|. The atomic sample is 
spin-polarized (mp = 1/2 or — 1 / 2 ) by properly detuned 
556 nm optical pumping in the presence of an 8.8 Gauss 
bias magnetic field. The inhomogeneity of our system 
is 0.23 at a temperature of 4.9(0.4) pK. To decrease the 
inhomogeneity, atoms in higher vibrational states of the 
optical lattice have to be blown out 0. We drop the 
intensity of the lattice laser to 53% level, at which the 
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FIG. 4: (Color online) Normalized collisional frequency shifts 
S co i as a function of the pulse area in Yb. The black dashed 
line (7 = 0.23), the red solid line (7 = 0.15), and the blue 
dashed-dotted line (7 = 0 . 10 ) are obtained from the analytical 
formula. The black squares and the red circles depict the 
experimental data at 7 = 0.23 and 7 = 0.15, respectively. The 
only s-wave collision case, S^ co i/U eg is plotted as the green 
dotted line. The inset shows the Rabi contrasts at 7 = 0.23 
(black dashed line) and 7 = 0.15 (red solid line). 

potential depth reaches 9 SE r . After 20 ms, the intensity 
of the lattice laser is increased to the original level (the 
intensity down-up method). The inhomogeneity is low¬ 
ered to 0.15, where the temperature is 3.9(0.3) fiK. The 
pulse area is controlled by the power of the clock laser at 
a fixed short pulse duration of T = 10 ms. 

Figure |4] shows the experimental and calculated nor¬ 
malized collisional frequency shifts as a function of the 
pulse area. The analytical solutions show good agree¬ 
ment with the experimental data. The only s-wave col¬ 
lision shift (the green dotted line) does not explain the 
measured shift and therefore we confirm that the p-wave 
collisions are dominant, which is consistent with the re¬ 
sult of Ramsey spectroscopy [13]. Although the mea¬ 
sured shift does not cross the zero shift mark due to the 
inhomogeneity limit of 0.15 in our experimental setup, 
the calculation proves that over-7r pulse cancellation is 
possible at an inhomogeneity level of less than 0.15, as 
depicted by the blue dashed-dotted line shown in Fig. 0J 
The lower inhomogeneity is attainable via an intensity 
down-up operation of the optical lattice, as described in 
our experimental description. 

The inset of Fig. 0] shows contrast of the Rabi oscilla¬ 
tion as a function of pulse area. The measured Rabi oscil¬ 
lation can be fitted by (A, 7) ~ (1 — cos Acos(^A))/2 
in the perturbative regime. We can extract the inhomo¬ 
geneity of the experiment. In the inset of Fig. 0J the 
maximum excited population of the over-7r pulse regime 
is smaller than that of the 7 r pulse. The S/N ratio of the 
excited population is proportional to y/(1 — p)/p where 
p is the excitation fraction [26|. Therefore it decreases 
due to the reduced excited population. 

The frequency noise of the clock laser and lock-servo 
error cause the population deviation from the 90% lock 


point, therefore induce the unwanted lock point variation 
which is a one of shift sources. The over-7r overdriving 
scheme is more sensitive to the frequency noise of the 
clock laser with respect to the under-7r scheme. In the ex¬ 
periment, the additional uncertainty ui oc k from the ±6% 
fluctuation of the 90% lock point leads to the residual 
frequency shift change ui oc k of a 10“ 17 level [25]. 

We repeat the five independent measurements at the 
IAtt pulse area. The A = IAtt pulse induces a collisional 
frequency shift closest to zero at an inhomogeneity level 
of 0.15. The shifts and uncertainties of the experiment 
are given in Table HT1 Comparing u co i of the A = IAtt case 
with that of the A = 0.67T case, we show that the over- 
7 r pulse gives the system a low level of total uncertainty 
by decreasing the uncertainty from the atomic density 
fluctuation at the expense of a small ui oc k increment. A 


TABLE II: Uncertainty budget of the experiment where V eg — 
V gg = -0.47 Hz, U ee - v gg = 0.39(Ue, - V gg ), U eg - V gg = 
-0.35(U e , - V gg ) 


A 

7 

$col 

^col 

'U'lock 

^tot 

'U'tot 



(xlO -3 ) 

(xlO -3 ) 

(xlO -3 ) 

(xlO -3 ) 

(xlO -18 ) 

0.6tt 

0.15 

638 

115 

12 

116 

105 

IAtt 

0.15 

74 

26 

19 

26 

29 


small fluctuation below ±1% and a small inhomogeneity 
below 0.1 provide the total uncertainty u to t less than 5 x 
10 -18 when \V eg — V gg \ < 0.6 Hz. 

In conclusion, we analyzed the collisional frequency 
shift in an optical lattice clock interrogated by a Rabi 
pulse. For the perturbative regime in which most op¬ 
tical lattice clocks are operated, the analytical solution 
was also obtained. It provides convenient expressions 
to explore the collisional shift and can be used to ex¬ 
tract collision interaction energies in Rabi spectroscopy. 
Based on our analysis, we propose that an over-7r pulse 
combined with a small inhomogeneity enables the cancel¬ 
lation of the total collisional frequency shift. This shows 
the potential for an optical lattice clock with a 10 -18 
uncertainty level with Rabi spectroscopy. 
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